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ABSTRACT 

We present three-dimensional global numerical simulations of the cosmic-ray 
(CR) driven dynamo in barred galaxies. We study the evolution of the interstellar 
medium of the barred galaxy in the presence of non-axisymmetric component of 
the potential, i.e. the bar. The magnetohydrodynamical dynamo is driven by 
CRs, which are continuously supplied to the disk by supernova (SN) remnants. 
No magnetic field is present at the beginning of simulations but one-tenth of 
SN explosions is a source of a small-scale randomly oriented dipolar magnetic 
field. In all models we assume that 10% of 10 51 erg SN kinetic energy output is 
converted into CR energy. 

To compare our results directly with the observed properties of galaxies we 
construct realistic maps of polarized radio emission. The main result is that the 
CR-driven dynamo can amplify weak magnetic fields up to a few //G within a 
few Gyr in barred galaxies. The obtained e-folding time is equal to 300 Myr and 
the magnetic field reaches equipartition at time t ~ 4.0 Gyr. Initially, completely 
random magnetic field evolves into large-scale structures. An even (quadrupole- 
type) configuration of the magnetic field with respect to the galactic plane can 
be observed. Additionally, the modeled magnetic field configuration resembles 
maps of the polarized intensity observed in barred galaxies. Polarization vectors 
are distributed along the bar and between spiral arms. Moreover, the drift of 
magnetic arms with respect to the spiral pattern in the gas density distribution 
is observed during the whole simulation time. 

Subject headings: cosmic rays - galaxies: evolution - galaxies: magnetic fields - 
methods: numerical 
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INTRODUCTION 



Bars are ubiquitous and occur in all types of disk galaxies, from early to late Hub- 
ble types. In near-infrared images about 70% of all the nearby disk galaxies are barred 
( Menendez-Delmestre et al. 2007). Bars are astrophysically important not only because 



they are very common in disk galaxies, but they also can significantly affect the gas dis- 
tribution (Sellwood & Wilkinson 1993) as well as the magnetic field configuration. The 



first systematic observations of polarized radio emission from a sample of 20 barred galaxies 
show that their magnetic field topology is significantly more complicated than in the case of 



grand-design spiral galaxies (Beck et al. 2002). The strongest polarized radio emission was 



detected in NGC 1097 and NGC 1365 ( |Beck et al.|2005| >. The main magnetic field properties 
observed in the above barred galaxies can be summarized as follows. The polarized emission 
is the strongest in the central part of the galaxy, where the bar is present. In this region the 
polarized emission forms ridges coinciding with the dust lanes along the leading edges of the 
bar. The polarization vectors change quickly their pitch angles in the bar region whenever 
they are located upstream the dust lanes and this results in depolarization valley where the 
polarized emission almost vanishes. Near the shear shock areas the regions of vanishing 
polarized intensity are also observed. In the outer disk magnetic field vectors form spiral 
pattern with the maxima of emission along spiral gaseous arms and in interarm regions. 
The average total (regular and turbulent) magnetic field strength for this sample of barred 
galaxies is 10 ± 3 fiG. On the other hand, the average regular magnetic field, calculated 
from polarized radio emission, is equal to 2.5 ± 0.8 fiG. The strongest total magnetic field 
is detected in the central starforming regions (about 60 /zG in NGC 1097) and in the radio 

ridges along the galaxies bars (20 /xG 30 fiG in NGC 1365). In the spiral arms of barred 

galaxies the total magnetic field is about 20 /iG, while regular 4 /iG. The spiral shape of the 
magnetic field, large pitch angles, and the observed strengths of magnetic fields indicate that 
the galactic dynamo works in those galaxies. 



The original concept of the fast CR driven dynamo was proposed by Parker (1992). Sev- 



eral researchers have approached this problem numerically, e.g., Hanasz et al. (2004, 2009a); 



Otmianowska-Mazur et al. (2009); Siejkowski et al. (2010). The CR-driven dynamo involves 



the following processes. CRs are continuously supplied to the galactic disk due to super- 
nova (SN) explosions. The galactic disk stratified by gravity is unstable due to buoyancy 
of the magnetic field and CRs. Buoyancy effects induce the formation of magnetic loops in 
the frozen-in, predominantly horizontal magnetic fields. The rotation of the interstellar gas 
causes that magnetic field loops are twisted by the Coriolis force. Next, due to the fast mag- 
netic reconnection small-scale magnetic loops merge to form the large-scale radial magnetic 
field component. The newly created magnetic field component is stretched by differential 
rotation, what results in amplification of the large-scale toroidal magnetic field component. 
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The combined action of the above effects is sufficient to trigger the exponential growth of 

the large-scale magnetic field with timescales 140 250 Myr, which are comparable to the 

galactic rotation period. 

Three-dimensional (3D) MHD numerical simulations in the shearing-box approximation 
have shown that the CR-driven dynamo can exponentially amplify weak magnetic fields 
up to few iiG within few Gyr in spiral galaxies (Hanasz et al. 2004 2009a) as well as in 
irregular galaxies ( Siejkowski et al.|2010 ). What is more, some of the observed magnetic fields 
properties such as extended halo structures of the edge-on galaxies or the so called X-shaped 
structures ( Soida|20"05 ), can be explained using the CR-driven dynamo ( Otmianowska-Mazur 



et al. 2009). The first complete global-scale 3D numerical model of the CR-driven dynamo 



in an axisymmetric galaxy has been demonstrated recently by Hanasz et al. (2009a). These 
simulations have given very interesting results and have shown that the CR-driven dynamo is 
one of the most promising process responsible for amplification and maintenance of galactic 
magnetic fields. 



2. THE INITIAL CONDITIONS AND INPUT PARAMETERS 



We investigated the evolution of the barred galaxy using the magnetized fluid approxi- 



mation governed by the isothermal non-ideal MHD equations (Landau & Lifshitz 1984). To 



make the set of above equations complete we include the CR transport. Following |Schlick- 



eiser & Lerche (1985), the propagation of CR component in the interstellar medium (ISM) 
is described by the diffusion-advection equation 

dc 

" ~ V(KVe cr )- Pcr (V -v) + CR source , (1) 



dt 



+ V(e cr v) 



where e cr is the CR density, p cr = {j cr — l)e cr is the CR pressure, K is the diffusion tensor, 
v is the velocity and CR source is the source term for cosmic ray energy. We assume that 
CR energy is added to the system by SNe explosions and that 10% of 10 51 erg SNe kinetic 
energy output is converted into CR energy, while the thermal energy is neglected. The 
adiabatic index 7 cr for the CR fluid is set to be 14/9. Additionally, CRs are weightless, 
thus they only contribute to the total pressure (not to the total mass) and are included in 



the gas motion equation as Vp cr (Berezinski et al.| 1990 ). Following Ryu et al. (2003), the 



anisotropic diffusion of the CR gas is described by diffusion tensor K as 

Kij = K ± Sij + (K\\ - K ± )n inj , (2) 

where K» and K± are parallel and perpendicular (with respect to the local magnetic field 
direction) CR diffusion coefficients and = Bi/B are components of the unit vectors tangent 
to magnetic field lines. 
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All numerical simulations were performed with the aid of the Godunov code (Kowal et al. 



2009). We numerically investigated the CR-driven dynamo model in 3D in a computational 



box of the size 30kpc x 30kpc x 8kpc, with a spatial resolution of 320 x 320 x 80 grid zones 
in the x, y, and z directions, respectively. The speed of sound is set to c s = 5kms _1 and 
the initial gas density at the galactic center po is equal to 1.0 at Hem -3 . Following several 



detailed reviews of the theory of CR diffusion (e.g. Strong et al. 2007), the values of the CR 
diffusion coefficients assumed in the simulations are Kn = 3 x 10 cm s" - 1 = lOOkp^Gyr" 1 



and K\ 



3 x 10 26 cm 2 s _1 



1 kpc Gyr 



-i 



whi 



our previous work (| Kulesza-Zydzik et al.||2009[ ) and it is set to be r\ 
0.1 kpc 2 Gyr" 1 



The resistivity coefficient is the same as in 

3 x lO^cmV 1 = 



e P — Per /p gas is constant and equal to 1 in the initial condition. At 



the beginning of the calculation the magnetic field is not present. Following Hanasz et al. 



(2009a), the magnetic field is added to the galactic disk through randomly distributed SN 
explosions in the period of time 0.1 Gyr 1.1 Gyr. During this time interval weak 10~ 5 fiG 



and dipolar magnetic field is supplied in 10% of SN remnants (Hanasz et al. 2009a). In 



the presented simulations the probability of a single SN event is proportional to the local 
gas density. After t = 1.1 Gyr the dipolar magnetic field is not longer injected because, 
due to the dynamo action, its contribution starts to be insignificant. We apply the outflow 
boundary conditions on external domain boundaries. 

Our model of the barred galaxy consists of four components: the large and massive 
halo, the central bulge, rotating disk of stars, and the bar. They are represented by different 
analytical gravitational potentials: the halo and the bulge components are described by two 
Plummer spheres, stellar disk is represented by the isochrone gravitational potential, while 



the bar is defined by the prolate spheroid (Binney & Tremaine 1987). The bar component 
is introduced into the galaxy gradually in time, until it reaches its final mass Mf, ar (from 
t = 0.1 Gyr to t = 0.4 Gyr). In order to conserve the total mass of the modeled galaxy we 
reduce the bulge mass, so we have Mt, ar (t) + M^it) = const during the calculations. The 
bar rotates with constant angular speed flbar- All quantities which characterize the model 
of barred galaxy are summarized in Table [TJ 



3. Magnetic field evolution 

In Figure [T] the toroidal magnetic field component in horizontal and vertical slices is 
plotted. Red color represents the regions with the positive toroidal magnetic field, blue with 
negative, while unmagnetized regions of the volume are white. The random magnetic field 
component dominates initially (up to about t ~ 1.0 Gyr), as it originates from randomly 
oriented magnetic dipoles. The ordered magnetic field is visible in the inner part of the 
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Parameter 


Name 


Value 


Units 


M d 


Disc mass 


4.0 x 10 10 


M e 


a d 


Length scale of the disk 


0.6 


kpc 


M b 


Bulge mass 


1.5 x 10 10 


M 


a b 


Length scale of the bulge 


5.0 


kpc 


M h 


Halo mass 


1.2 x 10 11 


M 


a h 


Length scale of the halo 


15.0 


kpc 


M bar 


Bar mass 


1.5 x 10 10 


M 


a bar 


Length of bar major axis 


6.0 


kpc 


bbar 


Length of bar minor axis 


3.0 


kpc 


Cbar 


Length of bar vertical axis 


2.5 


kpc 


^bar 


Bar angular velocity 


30.0 


Gyr- 1 


RgB 


Galaxy radius 


13.5 


kpc 


CR 


Corotation radius 


6.0 


kpc 


IILR 


Inner Inner Lindblad Resonance 


0.4 


kpc 


OILR 


Outer Inner Lindblad Resonance 


3.0 


kpc 


OLR 


Outer Lindblad Resonance 


8.5 


kpc 



Table 1: Adopted Parameters for the Barred Galaxy Model 
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Fig. 1. — Distribution of the toroidal magnetic field in vertical and horizontal slices through 
the disk center for selected time steps. Red color represents the regions with the positive 
toroidal magnetic field, blue with negative, while unmagnetized regions of the volume are 
white. To enhance weaker structures of the magnetic field in the outer galactic disk (e.g., 
magnetic arms), the color scale in magnetic field maps is saturated. 




Fig. 2. — Face-on and edge-on polarization maps for selected time steps. Polarized intensity 
(contours) and polarization angles (dashes) are superimposed onto column density plots. 
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galaxy, where it follows the gas distribution, namely the bar and dust lanes. At later epoch 
(t = 2.0 Gyr) the toroidal magnetic field component forms well-defined magnetic arms which 
can be observed till the end of the simulation. 

The positive toroidal magnetic field component gradually reaches higher values, both in 
the bar and magnetic arms, and occupies larger volume of space in the galactic disk and halo. 
On the other hand, the gas motions in the bar region generate shocks which enhance the 
magnetic field in this area in addition to the ongoing dynamo process. Indeed, the regular 
magnetic field in magnetic arms at time t = 5.5 Gyr is equal to 7.3 fiG, while in the bar 
region the regular magnetic field reaches 50.1 /zG. 

Reversals of the magnetic field can be observed at time t = 2.0 Gyr (Figure [T]) in the 
very inner part of the bar, where the negative toroidal magnetic field component is present. 
However, these reversals disappear almost completely in subsequent epochs. Additionally, 
due to the influence of the nonaxisymmetric gravitational potential, magnetic reversals are 
apparent between the bar and magnetic arms at time t = 5.5 Gyr in Figure [TJ 

At the beginning of the calculation the randomly distributed toroidal magnetic field 
dominates what is apparent in vertical slices. Next, at time t = 1.0 Gyr (Figure [T]), the 
odd (dipole-type) configuration of the magnetic field with respect to the galactic plane can 
be observed. However, this configuration is transient and after t = 2.0 Gyr in Figure [T] 
an even (quadruple-like) symmetry of the magnetic field dominates. The quadrupole-like 



symmetry was also obtained by Hanasz et al. (2009a) who studied the CR driven dynamo 



in an axisymmetric galaxies. Moreover, small reversals appear during the whole simulation 
time. 

In Figure [2] we present the magnetic field evolution in synthetic polarization maps for the 
same epochs as in Figure [T] The polarization maps show the distribution of the polarization 
angle and polarized intensity superimposed onto the column density. All face-on and edge- 
on polarization maps have been smoothed down to the resolution 40". The dark shades 
represent regions of the highest density. The magnetic field maxima correspond to the 
gas density enhancement, where the SN explosions are located, what can easily be seen at 
t = 1.0 Gyr (Figure [2]). At this time step, the magnetic field is present in the gaseous arms as 
well as in the central part of the galaxy. In the shock regions in the bar the rapid change of 
the magnetic field direction is apparent. Moreover, no regular magnetic field is observed in 
the interarm regions, however magnetic arms start to detach from gaseous spirals. At time 
t = 2.0 Gyr the magnetic spiral is clearly visible between the bar and gaseous arms. The drift 
of magnetic arms is persistent and takes place during the whole simulation time. Its shape 
changes slightly as the simulation proceeds. For instance, at t — 2.0 Gyr the magnetic spiral 
is well defined in the interam region, while at t — 5.5 Gyr it is less visible because it connects 
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with magnetic structures apparent in the bar region. The drift of magnetic arms into the 



interarm area was described in a number of papers e.g., (Kulesza-Zydzik et al. 



2009 



2010) 



where the authors concluded that this behavior is caused by difference in the angular velocity 
of magnetic arms and the gaseous spirals. Namely, the magnetic arms do not corotate with 
gaseous spiral structure, but they follow the general gas motion in the disk, which has a 



slightly lower angular velocity. However, in Kulesza-Zydzik et al. (2009, 2010) no dynamo 



action was included and the drift of magnetic arms into the interam region is not observed 
during the whole simulation time but only in the short period of calculation. 

In the edge-on maps in Figure [2j the polarized synchrotron emission extends far outside 
the galactic plane. Near the disk plane the magnetic field is mainly parallel do the disk, while 
in the halo vertical magnetic field components can also be seen. The strong vertical field 
in the halo is probably transported by mass outflow from the galactic disk. The averaged 
outflow rate is equal to l.lM yr _1 . The most extended structures are apparent at time 
t = 5.5 Gyr. This configuration of the magnetic field vectors bears some resemblance to the 



extended magnetic halo structures of the edge-on galaxies (Krause 2009). 



The CR-driven dynamo action in barred galaxies causes that the total magnetic field (left 
panel) and azimuthal flux (right panel) plotted in Figure |3]grow approximately exponentially. 
The growth of the total magnetic field lasts until the saturation level is reached at time 
t = 4.5 Gyr. The growth of the regular magnetic field is identified with the amplification of 
the azimuthal magnetic flux. The regular magnetic field grows up on an average timescale 
(e-folding time) equal to 300 Myr. After time t = 1.32 Gyr, the growth of the magnetic flux 
is clearly exponential and lasts until the equipartition is reached. 



4. Conclusion 

This work demonstrates for the first time the action of the CR driven dynamo in barred 
galaxies. Many observational futures of the magnetic field in barred galaxies have been 
reproduced. 

• The polarized radio emission is strongest in the inner part of the bar and in radio 
ridges that approximately follow the dust lanes indicated by the enhancement of the 
gas density. 

• In the outer part of the disk magnetic vectors form a spiral pattern with maxima of 
polarized intensity along spiral gaseous arms and in interarm regions. The drift of 
magnetic arms into the interarm area is observed during the whole simulation time. 
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• The synthetic edge-on radio maps of polarized emission resemble magnetic structures 
observed in edge-on barred galaxies. 

• The obtained average strength of the magnetic field and the maximum value of the total 
magnetic field in magnetic arms as well as in the bar are consistent with observational 
values. 

• The large-scale magnetic field has a quadrupole-like symmetry with respect to the 
galactic plane. 



The CR-driven dynamo is a powerful mechanism for amplifying magnetic fields in barred 
galaxies. During the whole simulation time we see the contribution driven by CRs driving 
turbulence, apparent through the buoyant motions. We incorporate the sub-grid physics 
through the resistivity of ISM. This amount of resistivity 77 = 3 x 10 25 cm 2 s _1 is exactly 
what is needed to ensure that magnetic energy in small-scale magnetic fields and in large- 
scale magnetic fields are equal in the CR-driven dynamo model (Hanasz et al. 2009b). The 
corresponding diffusive timescale is the magnetic field-ordering timescale. The timescale 
tgrow ~ \n(B / B )l / v tells us how much the random magnetic field on scales (I and v) would 
grow in the absence of resistivity. But since, due to the resistivity, the conversion of small- 
scale magnetic fields to large scales is efficient, we observe the growth of the large-scale 
magnetic fields to the equipartition values (from the initial B ), exactly in the period of 3-4 
Gyr. 



This work was supported by the Polish Ministry of Science and Higher Education 
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Fig. 3. — Evolution of the total magnetic energy Eb (left panel) and the mean azimuthal 
flux Bfj, (right panel) for the modelled barred galaxy. Both quantities are normalized with 
respect to the equipartition value. 



